Direct observation of a Fermi surface and superconducting gap in LuNi 2 B 2 C 
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We measured the Fermi surface (FS), band dispersion and superconducting gap in LUNI2B2C 
using Angle Resolved Photoemission Spectroscopy. Experimental data were compared with the 
tight-binding version of the Linear Muffin-Tin Orbital (LMTO) method and Linearized Augmented 
Plane- Wave (LAPW) calculations. We found reasonable agreement between the two calculations 
and experimental data. The measured FS exhibits large parallel regions with a nesting vector that 
agrees with a previous positron annihilation study and calculations of the generalized susceptibility. 
The measured dispersion curves also agree reasonably well with the TB-LMTO calculations, albeit 
with some differences in the strength of the hybridization. In addition, the spectrum in the super- 
conducting state revealed a 2meV superconducting gap. The data also clearly shows the presence 
of a coherent peak above the chemical potential, \i that originates from thermally excited electrons 
above the energy of 2A. This feature was not previously observed in the Lu-based material. 

PACS numbers: 74.70.Dd, 71.18.4-y, 71.20.-b, 71.27.+a 



I. INTRODUCTION 

Rare earth nickel borocarbides RNi 2 B 2 C ( R - rare 
earth) constitute an interesting class of material d 1 * 2 * 3 * 4 * 5 !, 
in which there is a competition and coexistence between 
superconductivity and magnetism. Amongst these com- 
pounds, nonmagnetic LuNi 2 B 2 C has the highest super- 
conducting critical temperature of 16.6 KM The boro- 
carbides exhibit a peculiar anisotropy of the supercon- 
ducting gap, the character of which is still under de- 
bate. It is believed that the gap is highly anisotropic 
in the t wo non-magn etic compounds LuNi 2 B 2 C and 
YNi 2 B 2 C P' 1 8 1 9 1 10 1 11 1 12 1 . Its symmetry was proposed to 
be s + g®, which is consistent with certain experimen- 
tal resultsP but an anisotropic s-wave symmetry has 
also been considered^. Other experimental data indi- 
cate that the gap in YNi 2 B 2 C has point nodes along the 
(100) and (010) directions^. LuNi 2 B 2 C crystallizes in 
a body-centered tetragonal structure with lattice param- 
eters a = 3.4639 A, and c = 10.6313 A 1 ^. Its crystal 
structure consists of Lu-C layers with Ni 2 B 2 sheets in 
between. Previously calculations reveal that LuNi 2 B 2 C 
is characterised by a large density of states (DOS) at 
the Fermi energy (Ep) orig inating mainly, but not ex- 
clusively, from Ni d electron s' 14 * 15 * 16 !. Another interesting 
feature is a flat band along the L-X direction just above 
E F . The Fermi surface (FS) topo graphy of LuNi 2 B 2 C 
was studied by ab-initio calculation s 17 * 18 * ^ . Band struc- 
ture calculations^- 7 -* revealed a pronounced maximum in 
the generalized electronic susceptibility at (~0.6a*, 0, 
0), where a* = 2ir/a and most likely arising from large 
nested regions of the FS. Moreover, phonon softening 
was observed in LuNi 2 B 2 C by means of inelastic neutron 
scattering for a range of wave vectors around (0.5a*, 0, 



0)P^l Interestingly enough the magnetic ordering, which 
was found in RNi 2 B 2 C compounds with magnetic atoms 
R = Er, Ho, Tb and Gd manifest similar modulation 
vector usually close to (0.55a*,0,0pl The first experi- 
mental studies of the LuNi 2 B 2 C (RNi 2 B 2 C) Fermi sur- 
face were performed by means of two-dimensional angu- 
lar correlation of electron-positron annihilation radiation 
(2D-ACAR) and the data were compared to the Lin- 
ear Muffin-Tin Orbital (LMTO), local density approx- 
imation (LDA) calculations^. Nested parts of the FS 
were found with a nesting vector corresponding to both 
the phonon softening and the magnetic modulation vec- 
tors. The fraction of the FS participating in nesting was 
determined to be 4.4 ± 0.59cp^*. That study was how- 
ever limited only to a rough "calipering" of the Fermi 
surface. Knowledge of the experimental band structure, 
Fermi surface and quasiparticle properties is deemed es- 
sential to understand the interplay of the various inter- 
actions in these materials, as it may shed new light on 
other phenomena such as anisotropic superconductivity, 
the role of phonon softening and the relationship be- 
tween the superconductivity and magnetic ordering in 
the borocarbides. It is also a pre-requisit for direct deter- 
mination of the alleged anisotropy of the superconducting 
gap in these materials. In this report we present angle 
resolved photoclcctron spectroscopy (ARPES) measure- 
ments of the band dispersion, Fermi surface and super- 
condcuting gap in the borocarbide with the highest T c , 
LuNi 2 B 2 C. The experimental results were compared with 
the tight-bounding LMTO (TB-LMTO) method and the 
full potential LAPW calculations. We found reason- 
able agreement with theory. The most significant dif- 
ference between the calculations and experimental data 
is the strength of the hybridization. We also deter- 
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FIG. 1: (Color online) Comparison between the Fermi surface maps measured by ARPES and the linear muffin-tin orbital 
(TB-LMTO) calculation, (a) Sketch of the first Brillouin zone for LUN12B2C. (b-d) ARPES mapping at the chemical potential 
for incident photon energies of 128.13 eV, 119.44 eV, and 102.98 eV, respectively, (e-g) The Fermi surface maps obtained by 
TB-LMTO calculations for constant k z values equal to (e) 0.2, (f) 0.15 and (g) 0.8 expressed in the units of T - Z distance. 



mined the superconducting gap to be 2.58 meV (extrap- 
olated for T=0), in good agreement with the gap ex- 
pected from the superconducting transition temperature 
(2A/k B T c = 2.78). 



II. EXPERIMENTAL 

LUNI2B2C single crystals were grown at Ames Labora- 
tory by means of a high-temperature flux techniqu e! 3 * 21 ! 
The plate-like crystals were cleaved in situ at pressures 
better than 3 x 10 -11 Tr to reveal and maintain fresh a-b 
surfaces. The Fermi surface and band structure mapping 
were performed at the 7.0.1 beamline at the Advanced 
Light Source, using a Scienta R4000 analyzer. The en- 
ergy and angle resolution were set at ~ 30 meV and ~ 
0.5 deg, respectively. The energy gap was measured with 
a Scienta 2002 analyser and He-I photon source (hv — 
21.2 eV), in which the overall energy resolution was set 
at 2 meV. The normal state data were measured at the 
Synchrotron Radiation Center using the PGM beamline 
and Scienta 2002 endstation, with the energy and angu- 
lar resolution set at ~ 13 meV and 0.25 deg, respec- 
tively. Tight-binding linear muffin-tin orbital calcula- 
tions were performed by the TB-LMTO program, ver- 
sion 47^21, and the Full-Potential Linearized Augmented 
Plane- Wave (LAPW) calculations were performed using 
the Wien2k package^. 



III. RESULTS AND DISCUSSION 

Band structure and semi-planar Fermi Surface cuts 
were determined for incident photon energies 128.13 eV, 
119.44 eV and 102.98 eV (Fig. lb-d), where the T point 
in the Brillouin zone (Fig. la) corresponds to normal 
emission of electrons along the (001) direction. A sin r/r 
correction term (r is the distance from the T point) was 
used to account for mapping of the momentum space onto 
the angular distribution of photoclcctrons. The ARPES 
process in 3D materials leaves some ambiguity as to the 
k z component of the momentum (perpendicular to the 
sample surface), because it is not conserved in the pho- 
toemission process. This is due to jump of the potential 
at the sample surface. From conservation of energy and 
remaining components of the momentum one can cal- 
culate the relative changes of the k z for various photon 
energies. To obtain the offset one needs to seek guid- 
ance from the band structure calculations and identify 
the high symmetry points in the dataP3 This allows es- 
timation of the k z offset. The change of the wave vector 
component k z (parallel to the c axis) between scans in 
Fig. lb (^=128.13) and lc (/w=119.44) was calculated 
from momentum and energy conservation to be 0.33 of 
the T-Z distance. Similarly the change of the wave vec- 
tor component k z between scans in Fig. lc (hv—l\%A^) 
and Id (hv=l02M) was 0.66 of the T-Z. The calculated 
Fermi surfaces were obtained for constant k z values by 
means of the TB-LMTO method and are shown in Fig. 
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FIG. 2: Dispersion of the conduction bands obtained with a photon energy of 102.98 eV (same as Fig. Id), compared with the 
TB-LMTO calculation, a) Fermi surface map with momentum cuts indicated by the solid lines, in which cut (c) and (i) pass 
through the F and X points, respectively, b) Fermi surface contours obtained by the TB-LMTO calculation for the value of k z 
corresponding to the data in panel (a). Panels (c-i): measured band dispersion along the cuts indicated in (a). Panels (j-p): 
calculated band dispersion along the cuts marked in panel (b). 



1 panels e-g. We estimated the values of the inner po- 
tential, Vo = 9.4 eV and the work function <f) — 4.6 eV, 
by comparing the high symmetry points between the cal- 
culated Fermi surfaces and the experimental data. This 
allowed us to determine the offset of the photon energy 
that corresponds to k z = 0. 

The Fermi surface maps for the incident photon ener- 
gies of 128.13 eV and 119.44 eV reveal large parallel parts 
of the FS with essentially the same nesting vector (spac- 
ing between the linear sections): k n = (0.59 ± 0.04) a* 
for Fig. lb and k n = (0.58 ± 0.04) a* for Fig. lc. k z 
is expressed in the units of the T-Z distance, where the 
r point corresponds to k z — 0. Although the full three 
dimensional FS was not determined in great detail, a sim- 
ilar nesting vector was found for different k z values which 
indicates that the FS likely has considerable nesting prop- 
erties for a wide range of k z values. The constancy of the 



value of the nesting vector between A; z =0.15 and -0.2 is 
also consistent with results of calculations. The spacing 
between the parallel segments of the Fermi surface pre- 
dicted by TB-LMTO calculation is between 0.54a* and 
0.55a*, the LAPW calculation results (not shown) are 
0.50a* and 0.57a*, respectively. The detected k n is very 
close to the theoretically predicted value obtained from 
the generalized susceptibility^. Our results also agree 
reasonably well with the nesting vector previously deter- 
mined via 2D-ACARP. 

The Fermi surface map obtained at 102.98 eV very 
closely resembles the calculated Fermi surface for k z = 
0.8a*. The overall shapes of the measured and calculated 
Fermi surface sheets (Fig. Id and lg) are very similar, 
however there is one significant difference. In the calcu- 
lations the four oval parts of the Fermi surface centered 
about T-Z are well separated in momentum space (Fig. 
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lg), while the data reveals that they actually are con- 
nected at the edges (Fig. Id). A lack of separation in 
the experimental data may indicate that the hybridiza- 
tion gap is overestimated in the calculations. These oval 
parts arise from the intersection of the electron and hole- 
like bands. Interestingly enough at the edges along the 
diagonal directions (e. g. 110) the bottom of the electron 
band and the top of the hole band appear to be pinned 
at the chemical potential, resulting in a characteristic 
"flower" shape. 

In Fig. 2 we plot the band dispersion data along a few 
selected cuts in momentum space obtained at an incident 
photon energy of 102.98 eV (Figs. 2c-i), along with a cal- 
culated (TB-LMTO method) band dispersion for k z ~ 
0.8a* (Figs. 2j-p). The agreement between the mea- 
sured and calculated band dispersion is rather good, es- 
pecially in the proximity of the chemical potential. In the 
corresponding TB-LMTO calculations (Figs. 2j-p), the 
same overall features are well reproduced, which shows 
the validity of the calculation in this material to a certain 
extent. This agreement also validates the assignment of 
k z values to the cuts measured at various photon ener- 
gies, which is very important when studying 3D materials 
with ARPES. The most significant difference is the hy- 
bridization gap, which is quite large in the calculations 
but its signatures are for the most part absent in the 
measured data. For example in Fig. 2 panels (k) and 
(1) the high and low energy branches form hybridization 
gap of about 200 meV at E = -0.3 eV, while in the corre- 
sponding measured data (panels (d) and (e)) the bands 
appear to disperse without a signature of the hybridiza- 
tion gap. One should consider if the observed features 
have any relation to the superconducting gap asymmetry 
and observed phonon softening in LuNi 2 B 2 C. It should 
be noted that band structure calculations^^ show a flat 
band lies very close to, but slightly above, the Fermi level. 
This feature was unfortunately not observed in our data 
due to the Fermi function cut-off. However, a higher 
DOS near the Fermi level would explain the large num- 
ber of scattered electrons observed with k vectors along 
(110) and phonon softening for the discussed wave vec- 
tors. Consequently this may lead to an anisotropy of the 
superconducting order parameter. This is in agreement 
with the results proposing for YNi2B2C that the super- 
conducting gap is larger just at (110) and diminishes or 
even has nodes along the (100) and (010) directions^. 

Given the above concern, we measured the energy 
gap in LuNi 2 B 2 C by partial angle-integrated photoelec- 
tron spectroscopy and compared with the normal state 
Fermi surface. The opening of the superconducting gap 
is clearly shown in Fig. 3. In order to determine the 
magnitude of the gap, the Dynes function^ was fitted to 
the symmetrized^ spectrum (Fig. 3b). The fitted func- 
tion yields the gap value of A = 1.5 meV for the sample 
at T = 11 K with the T parameter equal to 0.05 meV. 

The striking feature in Fig. 3c is the pronounced peak 
above the chemical potential. This peak arises from ther- 
mal excitation of electrons above the 2 A. This points to 
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FIG. 3: (Color online) Superconducting gap of LUM2B2C a) 
measured at T = 11±1 K, compared with the normal state 
at T = 40 K. b) The Dynes function (solid red line) with the 
parameters A = 1.5 meV and T — 0.05 meV fitted to the 
symmetrised spectrum (solid black circles), c) enlarged por- 
tion of superconducting spectra from (a) close to the chemical 
potential. 



high DOS just above fi which is consistent with the idea 
that the fiat band along (110) direction, a large part of 
which is slightly over the Fermi energy playing an impor- 
tant role in this anisotropic superconductivity. Similar 
peaks were recently reported in Y based borocarbides^. 
According to BCS theory the energy gap value at zero 
temperature (Ao) is 2.45 meV for T c = 16 K supercon- 
ductor. Our A value of 1.5 meV obtained at T=11K cor- 
responds to (A )=2.6 meV, in excellent agreement with 
the BCS predictions. 



IV. CONCLUSIONS 

We have performed measurements of the Fermi surface, 
band dispersion and supercondcuting gap for highest T c 
rare earth nickel borocarbide superconductor LuNi 2 B 2 C. 
The experimental data were compared with two different 
density functional calculations. The overall agreement 
between theory and measurement is good. In the experi- 
ment, large parallel FS parts spaced with the vector k n — 
0.59a* have been found for two different incident photon 
energies, which is a confirmation of the previous theoreti- 
cal predictions^ and earlier experimental studies^. The 
calculated FS confirms the existence of large nested parts, 
with a nesting vector in good agreement with the ARPES 
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results presented here. The superconducting gap was 
measured and we also observed a coherent peak above 
the chemical potential. This peak arises due to electrons 
being thermally excited above the energy of 2A. 
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